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Abstract 

Entangled photon-pair sources based on spontaneous parametric processes are widely used in photonic quan- 
turn information experiments. In this paper, we clarify the relationship between average photon-pair number 
and the visibility of two-photon interference (TPI) using those entanglement sources. We consider sources that 

Oh' 

generate distinguishable and indistinguishable entangled photon pairs, assuming coincidence measurements that 
use threshold detectors. We present formulas for the TPI visibility of a polarization entanglement that take 
account of all the high-order multi-pair emission events. Moreover, we show that the formulas can be approx- 
l— ~~ 1 imated with simple functions of the average pair number when the photon collection efficiency is small. As a 
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result, we reveal that an indistinguishable entangled pair provides better visibility than a distinguishable one. 
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1 Introduction 

Entanglement is an important resource in quantum information technologies, such as quantum key distribution 
(QKD) and quantum computation [T]. An entangled state of photons are especially important for distribut- 
ing entanglement over long distance. The most common way of generating entangled photon pairs is to use 
spontaneous parametric processes, including spontaneous parametric downconversion (SPDC) [21 El HI El E] and 
spontaneous four- wave mixing (SFWM) [?J [HJ El EH E] ■ These processes have been used for generating entan- 
gled states in polarization, energy-time, time-bin, and spatial modes. These entangled photon-pair sources have 
been used for a number of photonic quantum information experiments, such as the violation of Bell's inequality 
[U El [71 E| , quantum teleportation [T^] , entanglement swapping [T2] , QKD [21 Q~5] , and quantum logic gates 
|16j . In addition, entanglement generation in the 1.5- /xm telecom band is now being intensively studied with a 
view to realizing quantum communication systems over optical fiber. 

In previous experiments, two-photon interference (TPI) is often used to confirm the generation of entangle- 
ment. TPI visibility is a good measure of entanglement quality. It is well known that TPI visibility is limited 
by "accidental coincidences" . Entangled photon-pair sources based on spontaneous parametric processes emit 
multi-pairs with finite probabilities, and these multi-pair events cause coincidences with uncorrelated photons, 
which are called accidental coincidences. We cannot suppress the degradation in visibility caused by accidental 
coincidences completely unless we use photon number resolving single photon detectors with perfect photon 
collection efficiencies. In fact, most previous experiments used "threshold detectors", which are single photon 
detectors without a photon number resolution ability, and the photon collection efficiencies were much smaller 
than unity. Therefore, it is essential to take account of accidental coincidences when evaluating entangled 
photon pair sources. 

Since accidental coincidences are caused by multi-pair events, it is obvious that the probability of observing 
an accidental coincidence depends on the number distribution of the entangled photon pairs. The work by 
Luo and Chan [17] and Lim et al. [18] suggests that the pair number distribution becomes thermal-like when 
the photons generated in multi-pair emission events are quantum-mechanically indistinguishable. On the other 
hand, it is known that the distribution becomes Poissonian when the photons are distinguishable [T9l ED] . Even 
when pairs are generated in a single spatial mode, the actual distribution changes from a thermal-like to a 
Poisson distribution depending on the conditions of the pump pulse width, the length of nonlinear interaction, 
and the coherence time of the pairs. In any case, the visibilities improve when the average pair number is reduced 
by using a weaker pump. Therefore, many previous experiments used a very weak pump so that the accidental 
coincidences caused by multi-pair emission became negligible. However, in certain quantum information systems, 
it is sometimes better to use a relatively large average pair number, accepting the relatively large number of 
accidental coincidences. For example, in an entanglement-based QKD system, we can optimize the average pair 
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number to maximize the secure key rate: i.e. we can often increase the secure key rate by using a relatively 
large average pair number while accepting an increased error rate. Since the error rate is directly linked to the 
TPI visibility, it is meaningful to know the relationship between the average pair number and the visibility. 

The purpose of this study is to clarify the relationship between the average pair number and the quality 
of the entangled photon pairs generated by spontaneous parametric processes. Although the TPI visibility of 
distinguishable entanglement sources is considered briefly in 19, 20J , its full derivation has not yet been reported. 
Moreover, the TPI visibility of an indistinguishable entanglement source has not yet been clarified. Therefore 
in this paper, we present the formulas for the TPI visibilities of distinguishable and indistinguishable entangled 
photon-pair sources. We assume the use of threshold detectors, since most of the entanglement generation 
experiments conducted so far have used this type of detector. In particular, when the collection efficiency 
of the photon is small (this assumption can be validated for most of the previous experimental conditions), 
the TPI visibilities are given by simple functions of the average pair number. Then, the density matrices 
and concurrences of distinguishable and indistinguishable entanglement sources are given as a function of the 
average pair number. We also present formulas modified for time-bin entanglement measurements. We believe 
that our results will be useful for designing and evaluating quantum information experiments using entangled 
photon-pair sources based on spontaneous parametric processes. 

The organization of this paper is as follows. In section we review the pair number distribution of indistin- 
guishable and distinguishable entangled photon pairs, which become thermal-like and Poissonian, respectively. 
In section [J2 the TPI visibilities of two types of entanglement sources are given as a function of average pair 
number and collection efficiency Section fj] shows that the density matrices of those sources can be theoretically 
obtained if the average pair number and the collection efficiencies are known. Here, we also calculate the theo- 
retical concurrence using the obtained density matrices. Section [5] considers correction terms that take account 
of the detector dark counts. Section [5] describes formulas that are modified for time-bin entangled photon pairs. 
Finally, the obtained results are summarized in section [7] 

2 Number distribution of entangled photon pairs generated by spon- 
taneous parametric processes 

In this section, we derive the pair number distributions when the generated multi-pairs are quantum-mechanically 
indistinguishable, and distinguishable. 
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2.1 Indistinguishable entangled photon pairs 

When photon pairs are generated in the same mode by a spontaneous parametric process, the photons obtained 
in a multi-emission event are indistinguishable. According to [21] . an interaction Hamiltonian for single-mode 
correlated photon pair generation by spontaneous parametric processes at the origin of the time evolution can 
be expressed as 

H mt = H%x{a\a\ - a s a,-), (1) 

where subscripts s and i denote the signal and idler modes, respectively, x m the above equation for SPDC and 
SFWM is expressed as 

Xspdc = cx {2) a p , (2) 
Xsfwm = c'x (3) a P ia P 2, (3) 

where a p , a p \ and a v i denote pump field amplitudes. Here, we assume that the pumps are strong and can be 
treated as classical oscillators. 

Entangled photon pairs are generated when two spontaneous parametric processes in different modes occur 
coherently. Here we assume polarization entangled photon pairs, which are generated by the following interaction 
Hamiltonian. 

H-i n tE = HintH + H in tV ■ (4) 

Here, H and V denote the horizontal and vertical polarization modes, respectively, and Hi n tH and Hi n tv are 
the interaction Hamiltonians for generating H and V polarization photon pairs given by 

Hi n tH = i^x{a\ H a\ H - a sH a iH ), (5) 
Hi n tv = iftx( a lv a lv - a sV a iV ). (6) 

Using Hi n tE, the time evolution of state is described with the following equation. 

|*)=exp|-t^t||0) (7) 
According to the disentangling theorem shown as Eq. (5.63) in [21], the exponential part can be rewritten as 

. Hi n t 



exp 
where 



j intE I \ = e Ta l H AH e ~3(al H a s H+a\ H a iH + l) e -Ta sH a iH e Ta\ v a\ v e -g(a\ v a sV +a\ v a iv + l) e -Ta BV a iv ^ 



r = tanh%t, 
9 = ln(coshx£). 



(9) 
(10) 
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With this equation, and discarding the term that becomes zero, Eq. ([7]) is expressed as 

I*) = ^exp(iJ 2 )|0) 

= ^{ 1 + ^ + ^2 2 + --- + ^ + -4l°>- (11) 

Here, H 2 is defined as 



H 2 =T(al H a\ H +al v al v ), (12) 

and C = coshx^- 

The first order term in Eq. pip gives an entangled photon pair state, which is calculated as 

= -^H 2 \0) = ^{|l,0).|l,0) < + |0 > l).|0,l) i } J (13) 

where |o, b) x denotes that in mode x, there are a photons with H polarization and b photons with V polarization. 
Similarly, the second order term, in which two-pair state is generated, is expressed as follows. 

|*w»kK.) = -^Hl\0) = ^{|2,0) fl |2,0) i + |l,l) s |l,l>i + 10,2)^0,2),} (14) 

This is the state that is described as an entangled two-photon polarization state in [52] . 
The xth order term can be expressed as 

pa: x j~*x 

Y&xindis) = £Q./Jx- fc) k) s \x - k, k)i = V%+ l-^\ip x ), (15) 

fc=0 

where a; is a positive integer and is a normalized quantum state of x-pair given by 

1 x 

Wx) = —j — — ^ \x - k, k) s \x - k, k)i. (16) 
V^+l^ 

Equation (|15|) shows that the probability of generating x pairs Pindis (x) is given by 

(x + l)T 2x _ (x + 1) tanh 2 * X t 
C 4 cosh 4 xt 

and the average photon-pair number [i is given by 



. . _ -r ±j± _ -r ±j uaiiii x L / 17 n 

*tndia\ x ) — — ,4 , 1 



M = 2J xP indis(x) = 2 sinh 2 xt- (18) 

31=0 

When the nonlinear medium is pumped by pulses, the average pair number is defined as a number per pulse. 
When it is pumped by a continuous light, the average pair number is defined as a number per temporal resolution 
of a measurement system, which is often determined by the temporal resolution of the photon detectors. 

Using a relationship cosh 2 xt — sinh 2 xt = 1, the probability of generating x pairs is expressed as a function 



of /x, which is 
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2.2 Distinguishable entangled photon pairs 

We then consider the situation in which a pump can generate entangled photon pairs in many different modes. 
This situation is typically observed when a continuous wave pump is used or the pump pulse width is much longer 
than the coherence time of the photon pairs [531 HI]- I n those cases, the generated pairs are distinguishable 
in the time domain. Also, even if the pump pulse width is much shorter than the photon pair coherent time, 
the group delay difference between pump pulse and photon pair in the nonlinear crystal can lead to temporal 
distinguishability. In those cases, the pairs generated in multi-emission events are mutually distinguishable. 
For example, when two pairs of entangled photons are generated through independent spontaneous parametric 
processes in temporal modes 1 and 2, the joint state is given by 

|* 2 rf») = ^{|l,0) a i|l,0)« + |0 > l).i|0,l)fl}®{|l,0)«a|l > 0)« + |0,l) i 2|0,l) i2 }. (20) 

In the following, we explain why the number distribution of distinguishable entangled photon pairs becomes 
Poissonian. 

Here we consider the temporal distinguishability that arises from the use of a broad pump pulse. We define 
an integer TV as follows. 

N = At P um P (21) 
At pp 

Here, At pump and At pp denote the pump pulse width and coherence time of the photon pairs, respectively. The 
physical meaning of N is the number of different temporal modes excited by the pump pulse. Let us consider a 
case where x pairs are created. If N is far larger than x, the probability that two or more pairs are in the same 
temporal mode is negligible. In this case, the x pairs are all distinguishable. This means that the probability 
of generating x pairs is simply given by multiplying of the single pair generation probability, which is given by 
Eq. (|17p with x = 1. The number of combinations for choosing the temporal modes of the x pairs in N modes 
is given by nC x . Therefore, the probability of generating x pairs is expressed as 

[2T 2 \ X [2T 2 N\ X 1 , N 

P v {N,x)~a. N C x [^j ^a^-^J -, (22) 

where a is a constant for normalization, and an approximation nC x ~ N x jx\ is used. By setting ( 2V c i ^ = M 
and a — e~ M , the probability of generating x distinguishable pairs is obtained as 

P p (»,x)~e-^. (23) 
x\ 

Thus, when N >> x, the pair number distribution becomes Poissonian. It is obvious that /i gives the average 
pair number. 
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Figure 1: (a) TPI visibility measurement setup, (b) typical TPI waveform. 

3 Two-photon interference measurement 

In this section, we consider the visibility of the TPI measurement. We assume polarization entangled photon-pair 
states whose first-order approximations are given by 



i 



^=(|l,0) s |l,0) l + |0,l) s |0,l> 4 ) 



(24) 



In the following, we analyze the way in which the multi-pair emission degrades the visibilities of TPI fringes for 
distinguishable and indistinguishable entangled photon pairs. In this analysis we discard all the other causes of 
visibility degradation such as a deviation of the polarization reference between the signal and idler measurements 
and the dark counts of the detectors. 

We assume the experimental setup shown in Fig. [T] (a). A polarization entangled photon-pair source 
generates a signal and an idler photon. Each photon is input into a rotatable polarizer that is followed by a 
threshold detector. When x photons are incident on a threshold detector with a collection efficiency of a, the 
detection probability is obtained as 1 — (1 — a) x . Here, the collection efficiency is a value that includes all the 
optical losses and the quantum efficiency of the detector. The detection signals from the detectors are input into 
a time interval analyzer (TIA) to measure the coincidence counts. We assume that the polarizer for the signal is 
fixed at H polarization, and then the coincidence rate is measured as a function of idler polarizer angle. Figure 
[T] (b) shows a typical TPI fringe for the (approximated) state given by Eq. (j2"4"|) . where the coincidence rate 
is illustrated as a function of idler polarizer angle from H polarization. When we set the idler polarizer angle 
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at H and ^polarizations, we observe the "peak" and "bottom" of the TPI fringe, respectively. We denote the 
coincidence rates at the peak and bottom as Rhh and Rhv , respectively. Then, the visibility of the waveform 
is obtained as 

v _ Rhh — Rhv 
Rhh + Rhv 

We also introduce a contrast parameter r defined as 

r = ^ (26) 

In the following, we theoretically calculate Rhh and Rhv to obtain the v and r values for the distinguishable 
and indistinguishable entangled photon pairs. 

3.1 Distinguishable entangled photon pairs 

Here, we consider the TPI visibility for the distinguishable entangled pairs whose pair number distribution is 
given by Eq. 



3.1.1 Coincidence rate at fringe peak 

We assume that both the signal and idler polarizers are set at H polarization to obtain coincidence rate Rhh ■ 
First, let us consider a case in which a single pair is created with a probability of P p (p, 1) given by Eq. 
(|23|) . The patterns of polarization for a photon pair are only both H or both V, which occurs with the same 
emergence probability of 1/2. When both photons have H polarization, a coincidence count is observed with a 
detection probability of a 2 . Therefore, the overall coincidence probability when one pair is created is given by 

Table 1: Single pair 



Pattern 


Emergence probability 


Detection probability 


HH 


1 

2 


a 2 


VV 


1 

2 






Then, we consider a case where two distinguishable pairs are created. The state is given by Eq. (|2"U|) . The two 
pairs are created with a probability of P p (fJ-, 2). Here we should note that the two pairs are distinguishable, so 
we denote them pair 1 and pair 2. The patterns of polarization, the emergence probabilities of the corresponding 
patterns, and the probabilities of detecting coincidences for the corresponding patterns using threshold detectors 
are summarized in Table 2. The table indicates that the total coincidence probability of two pairs is given by 
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Table 2: Two pairs 



Pattern 


Emergence probability 


Detection probability 


EE 
EE 


1 

4 




EE 


1 




VV 


4 


a 2 


vv 


1 




EE 


4 


a 2 


VV 

vv 


1 

4 






Then let us consider the ic-pair case, where a; is a positive integer. A distinguishable s-pair state is expressed 



as 



I*, 



fc + |0,l) sfe |0,l) ife }. 



(27) 



The total number of patterns is given by 2 X , and so the emergence probability of each pattern is given by 1/2 X . 
If y pairs are in the VV state (y: nonnegative integer), the number of patterns of different sets of EE and VV 
pairs is expressed as x C y — y \^L y y i with a detection probability of {1 — (1 — a) x ~ v } 2 . Thus, the coincidence 
probability in the setup when x pairs are generated is given by 



y=o 



2 X y\(x - y) 



■{i-(i-a) x -y} 2 



(28) 



Since the probability of observing x photon pairs is given by P x (fi), the overall coincidence rate is obtained 
with the following equation. 

oo 

fe P = ^/(^(f,^) (29) 

When a << 1, the result shown in Eq. (|2"9"|) can be approximated as follows. Using {1 — (1 — a) x ~ v } 2 ~ 
(x-y) 2 a 2 , 



OO ( X 



Rhh p 



EE 



i 



— (x — y) 2 ——e^ 

, ^ 2 X ylfx -y)\ 1 V > x\ 

x=0 (,y=0 " v y ' 



2 ( M , I* 



(30) 



The detailed derivation of the above equation is shown with Eqs. (j9"5f to (|100p in Appendix [Bl 
3.1.2 Coincidence rate at fringe bottom 

Here we assume that the signal and the idler polarizer are set at E and V polarizations, respectively, to obtain 
the coincidence rate Rhv- As with the "peak" case, when x pairs are created, the total number of patterns is 
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given by 2 X , and the emergence probability of each pattern is given by l/2 a: . If y pairs are in the VV state, the 
number of patterns of different sets of HH and VV pairs is expressed as x C y = y \ix—y\\ - Since the numbers of 
FUand HH pairs are y and x — y, respectively, the detection probability is given by {1 — (1 — a) x ~ v }{l — (1 — a) v }. 
Thus, the coincidence detection probability in the setup when x pairs are generated is given by 



= E h^vv {1 - (1 - a)x ~ V}{1 - (1 - a)v} (31) 

y=Q y ' V y) ' 

The overall coincidence rate is obtained with the following equation. 

OO 

Rhv p = ^2 9( x ) P p(V, x ) (32) 

x=0 

As with the peak coincidence rate, we can simplify the above equation when a << 1. Using {1 — (1 — 
a) x ~ v }{l — (1 — a) y } ~ (a; — y)ya 2 , Rhvp is reduced to 

- (33) 

In the transformation of the above equation, we used a similar procedure to that used to derive Eq. (|30|) with 
the relationship given by Eq. (jlOip . 

3.1.3 Single count rate 

If the idler detection information is not available to the signal side, the state of the signal photon can be 
expressed as a statistical mixture. When x pairs are generated, the state of the signal is expressed as 

^(g)(|l,0) sfc (l,0| sfc + |0,l) sfe (0,l| sfc ). (34) 
fc=i 

The state of the idler photon can also be expressed with a similar formula. If we assume that the polarizer for 
the signal is set at H polarization, the probability of detecting the above state, s(x), is given by 

s{x)=Y^ t/ Xl - {1 - (1 - a) g - y }. (35) 
Using s(#), the single count rate R sp is expressed as 

OO 

R sp = p p(p> x ) s ( x ) ( 36 ) 

x=Q 

With the approximation 1 — (1 — a) x ~ xa, the single count rate is reduced to 

Rsp -/xa. (37) 
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3.1.4 Visibility and contrast 

The visibility and the contrast with no approximation are obtained by plugging Eqs. 
following equations. 

Rhh p — R-HVp 
v p — - 

Rhh p + Rhv p 
_ Rhh p 



and 



into the 

(38) 
(39) 



Rhv p 

With 1 — (1 — a) x ~ xa, we can obtain the approximated visibility and contrast using Eqs. (l3^|) and (l33|) as 

1 



1+M' 
2 

1 + -. 

A 4 



(40) 
(41) 



Thus, we can ignore the loss dependence of v p and t p when a << 1. It is obvious from the derivation shown 
above that the above approximated equations take account of both the limited number of multi-pair emission 
events and all the high order multi-pair emission events. 

3.2 Indistinguishable entangled photon pairs 

Here we calculate the visibility of the TPI fringe obtained with indistinguishable entangled photon pairs, whose 
pair number distribution is given by Eq. (|19[) . 



3.2.1 Coincidence rate at fringe peak 

The detection probability when a single pair is created is the same as that for distinguishable pairs. However, 
the emergence probabilities of multi pairs are different from those for distinguishable pairs, since the number of 
patterns of polarization decreases because of the indistinguishability of the pairs. For example, let us consider 
the case when two pairs are created. There are only three patterns: both are HH, one is HH and the other is 
VV, both are VV. The patterns are summarized in Table 3. 

Table 3: Two pairs 



Pattern 


Emergence probability 


Detection probability 


Both HH 


l 

3 


U-(i-") 2 } 2 


HH and VV 


1 
3 


a 2 


Both VV 


1 
3 






In general, when an x-pair state as shown by Eq. (|16| is created, the total number of patterns is given 
by x + 1, and so the emergence probability of each pattern is given by l/(x + 1). If y pairs are in the VV 
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state, the number of patterns of different sets of HH and VV pairs is simply 1, with a detection probability of 
{1 — (1 — a) x ~ y } 2 . Thus, the coincidence probability in the setup when x pairs are generated is given by 

hnM = r-rri 1 - a - a T~ y } 2 ( 42 ) 

Since the probability of observing x photon pairs is given by Pi n dis(n,x) (Eq. (|19p). the overall coincidence 
rate is obtained with the following equation. 



RHHindis = ^ ^ findis {p^)Pjndis (Mi ^) (43) 
x=0 

When 1 — (1 — a) x ~ ra, the above equation can be simplified as with distinguishable pairs. 

Ri 

In the transformation of the above equation, we used Eqs. ([M]) . and (|102p . 



iHHindis ^ X! jX! ~~ y) 2 a 2 Pindis{Hi x ) j = ^ + yj ^ 



3.2.2 Coincidence rate at fringe bottom 

We now consider the coincidence rate at the bottom of the fringe, Rhv- With a similar calculation, the 
coincidence probability in the setup when x pairs are generated is given by 



QinUx) = J2 ^T^ 1 ~ 0- ~ - (1 - «)»} (45) 

y=o 

The overall coincidence rate at the fringe bottom is obtained with the following equation. 

oo 

RHVindis = ^ dindis (x)Pindis (M, x) (46) 

We can also derive a simpler expression for RHVindis when 1 — (1 — a) x ~ ia. 

00 f x 1 1 2 a 2 

RHVindis - ^ < ^ — —(a; - y)ya 2 P indis (n, x) \ = — — (47) 

a;=0 U=0 X J 

In the above equations, we used the relationship given by Eqs. (|94|) . (|96|) and (|103|) . 
3.2.3 Single count rate 

When indistinguishable x pairs are generated, the signal photon state can be expressed as the following statistical 
mixture. 

—^—^2\x-k,k) sk (x-k,k\ xk (48) 

X k=0 
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If the above state is detected with the polarizer set with H polarization, the detection probability Si n dis{x) is 
given by 

s indis (x) = ^—J2{l-(l-a) k }. (49) 
The single count rate R S indis is given by 

oo 

Rsindis — ^ Pjndis (/^; ^Q^irerf^a (*^) (50) 
x—0 

With the approximation 1 — (1 — a) x ~ ra, the above equation is simplified to 

Rsindis - ^ | ^d™^" X ^ Y^ka) = i/ia (51) 
x=0 \ k=o / 

Thus, when the collection efficiency is small, the single count rates of the indistinguishable and distinguishable 
pairs are the same. 

3.2.4 Visibility and contrast 

The visibility and contrast for an indistinguishable source are given by 



Rh Hindis RhV indis { KO\ 

Vindis = "5 —5 (OZ) 
i^HHindis + RHVindis 

Rhh indis /r-o\ 

Tindis = "TJ \P") 
^HVindis 

When 1 — (1 — a) x ~ xa 7 the above values are simply obtained as 

Vindis - 3^ (54) 
2 

Tindis — 2 H (55) 

M 

Now we evaluate the validity of the visibility obtained with an approximation of 1 — (1 — a) x ~ xa. Assuming 
/j = 0.1, Eqs. (fS"2"|) and (j5"4")l are plotted as a function of collection efficiency a. The result is shown in Fig. 
Here, we took account of up to 100 terms in the calculation of Eq. (|5"2"]) . Although the deviation becomes larger 
as the collection efficiency increases, the difference in the visibility is only 0.4% even with perfect collection 
efficiency. In experiments, it is very hard to obtain a collection efficiency exceeding 50% even in the short 
wavelength band. Moreover, in most of the 1.5-/im band experiments, the collection efficiencies are typically a 
few per cent, largely because the single photon detectors for this band are inefficient. Also, we note that the 
difference between Eqs. (fS"2"T) and ([54"]) becomes even smaller if we use a smaller /i. Similar results are obtained 
with the formula for the distinguishable pairs. Thus, Eqs. (|40j) and (|54|) provide very good approximations of 
the visibilities. 
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Figure 2: Visibility as a function of collection efficiency using Eqs. ([52"|) (dotted line) and (|54"|) (solid line). 
3.3 Comparison of distinguishable and indistinguishable pairs 

We compared the visibilities of distinguishable and indistinguishable entangled photon pairs using Eqs. and 
(Jill), where we assumed a = 0.1 and took account of up to 100 terms in the calculations of these equations. The 
result is shown in Fig. [31 This indicates that the indistinguishable pairs exhibit slightly better visibility than 
the distinguishable pairs. For example, at fi = 0.1, the visibilities for the indistinguishable and distinguishable 
pairs are 91.2% and 90.8%, respectively. We also note that we can obtain almost the same results if we use 
approximated equations (|4T)|) and ([54"|) . The result reflects the fact that the contrast of indistinguishable pairs 
Tindis is greater than that of distinguishable pairs t p by ~ 1, as we saw in the approximated equations of the 
contrasts given by Eqs. (|4"Tj) and In addition, Eqs. ([50)1 and (f4"4"]) suggest that the indistinguishable pairs 

yield a larger peak coincidence rate than the distinguishable pairs. These results mean that, in an entanglement- 
based QKD system, an indistinguishable entangled pair source gives a larger key rate with a smaller error rate 
than a distinguishable source, which yields a higher secure key rate. Thus, the presented analysis clearly 
confirms that an indistinguishable entangled photon pair source performs better than a distinguishable one in 
quantum information experiments. 
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Figure 3: Visibility as a function of average photon number per pulse. Solid line: distinguishable pairs (Eq. 
([55])), dotted line: indistinguishable pairs (Eq. (JHU)). 

4 Density matrix reconstruction using quantum state tomography 

The method for calculating the coincidence rate described in the previous section can be used to calculate 
the 16 different combinations of two-photon projection measurements shown in Table 4. In this table, \H) 
(\V)) denotes a state in which there is a photon in the H (V) polarization mode. Other states are defined as 
|+) = (\H) + \V))/y/2, \L) = (\H)+i\V))/V2, and \R) = (\H)-i\V))/V2. By using these coincidence rates, we 
can construct the density matrix of the two-photon state using quantum state tomography (QST) as reported 
in |25j . Here, again we assume that only the accidental coincidences caused by multi-photon emission events 
are the origin of the degradation in the entanglement quality, and there is no other cause for degradation. 
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Table 4. Assumed projected state. 



V 


Photon 1 


Photon 2 


1 


\H) 


\H) 


2 


\H) 


\V) 


3 


\V) 


\V) 


4 


\V) 


\H) 


5 


\R) 


\H) 


6 


\R) 


\V) 


7 


1+) 

1 1 / 


\V) 


8 


1+) 

1 1 / 


\H) 


9 


1+) 

1 1 / 


\R) 


10 


1+) 

1 1 / 


1+) 

1 1 / 


11 


\R) 


1+) 


12 


\H) 


1+) 


13 


\V) 


1+) 


14 


\V) 


\L) 


15 


\H) 


\L) 


16 


\R) 


\L) 



v = 1, 3, 10 and 16 correspond to positively correlated events. Therefore, it is easy to show that these 
events have the same coincidence rate as Rhh calculated in the previous section. Similarly, the coincidence 
rates for v =2, 4, which correspond to negatively correlated events, are given by Rhv- From the symmetry 
of the equations, the coincidence rates for other basis mismatch cases have the same value. Therefore, if we 
can calculate a coincidence rate for one of those combinations, we can obtain the rates for all v values and 
thus construct the density matrix. In the following, we consider the v = 12 case, in which the signal photon is 
projected to the \H) state and the idler photon is projected to the |+) state. 

4.1 Distinguishable pairs 

Assume that the projection measurement for the signal photon is with H V basis, and that for the idler is with 
diagonal basis. The single pair state is rewritten as 

= -^(\H)s\H) t + \V) s \Vh) 

= lQH) a \+) i + \V) B \+) i + \H) a \-) i -\V),\-) i ). (56) 
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When an x-pair state given by Eq. (|27p is generated, the coincidence probability is given by 



h(x) = 



x 1 
oa 



\k=Q 



2 X k\{x-k)\ 



{i -(i -ay} 



(57) 



The overall probability of coincidence is obtained as 

oo 

r v is then obtained using coincidence rates of HH, HV and H+ projection measurements. 

Rhh p 1^=1,3,10,16 
Rhvp v — 2, 4 
Rh+p v — others 

When 1 — (1 — a) k ~ ha, Eq. (|58|) can be simplified to 

2 ( V , 



(58) 



(59) 



(60) 



4 4, 

In the following, we calculate the density matrix assuming 1 — (1 — a) k ~ fca. Let us denote the coincidence 
rate for each value in Table 4 as r„. Using the results of the above equation as well as those of the previous 
section, r„ is obtained as follows. 



a 2 (f + ^J i/ = l > 3,10,16 

2 2 



(61) 



Using r v shown above and the procedure presented in |25j , we can obtain the following density matrix. 



( 2 +m 

4+4/x 






4+4 M 









\ 



4+4 M 







4+4^ / 



(62) 



V 2+2p 

Thus, the density matrix is obtained as a function of ji. The density matrix is graphically shown in Fig. H] for 
/i = 0.01 and 0.3. Using this density matrix, the concurrence can also be obtained as a simple function of //, by 
using the procedure described in [25j . The concurrence is 

2 - (i 



Conc n 



2(m + 1)' 



(63) 
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Figure 4: Graphical representation ol density matrix of distinguishable photon-pair states, (a) fi — 0.01, (b) 
fi = 0.3. 

4.2 Indistinguishable pairs 

Here we undertake a projection measurement on the indistinguishable pairs in the way described in the previous 
subsection. 

Let us assume that we have an x indistinguishable pair state, whose state is given by Eq. (|16p . In that joint 
state, \x — k, k)i can be rewritten as follows. 

1 



\x-k,k)i = - 7 =^==(at H r- k (4 v ) k \0) i 

1 

W2)^{x~k)\k^ " ' "' 



(aL+oJ_) B -*(aL-a{_)*|0) 



t \fci 



Here, we used the unitary transformations between HV and the diagonal basis given by 

a H = -^(a+ + al), 

a f v = -y=(a\. - ai). 
Using binominal theorem (a + b) k = ^ n=0 .kC n a n b k ~ n , Eq. is further transformed as 



(64) 

(65) 
(66) 



i 



£ . x . k C m {a\ + r{a\_) 



t \x — k- 



Y,-kC n {a\ + Y{-aU k - n )\Q)i (67) 



m^{ X -k)\k\ Lm=0 ) u=0 ) 

From the above equation, the coincidence probability of the x indistinguishable pair state is obtained as 



Hindis yfti OS.) 
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1 33—1 

— Y 



c — k k 

EE 



{(x-fc)!} 2 (fc!) 
2 x (x - k)\k\ ^ ^ m\(x - k - n)\ nUk - n)\ 

y ' m=On=Q y J \ J 



{1 - (1 - a) n+m }{\ - (1 - a) x ~ k } 



(68) 



Using hindis, we can obtain the overall coincidence rate for the \H) S and |+)^ bases with the following equation. 

oo 

(M. ") = E P indis (M; x)h lndis (p, a, x) (69) 



x=0 



r v is then determined using coincidence rates obtained with HH, HV and H+ projection measurements. 



RtlHindis v = 1, 3, 10, 16 
TV = \ RHVindis V = 2, 4 

Rn+mdis v = others 



(70) 



The density matrix, and consequently the concurrence, arc obtained using those rates as a function of fj, and a. 
When 1 — (1 — a) k ~ ka, Equation (p)9"|) is approximated as follows. 



3 M 2 



Ru+indisip, a) ~ a" 1 (^ _ ^ _ + ^ 
By using the approximated coincidence rate equations, r„ is now obtained as 

a 2 (4 + f) f= 1,3, 10, 16 



(71) 



(72) 



a 



{ 8 f i) 

With those r„ values, the density matrix is obtained as 

/ 



v = others 



Pindis 



1+fJ. 








2+M 


2+3/j 


4+6^ 













4+6p 








A* 





4+6^ 










1+H 


4+6p 


2+3p, 



(73) 



Using the above density matrix, we approximate the concurrence with the following equation. 

2 



Conctiji) 



(74) 



2 + 3/i 

Using Eqs. (f6"3"|) and (fT4")) . we plot the concurrences of distinguishable and indistinguishable entangled photon 
pairs as a function of /i in Fig. \E\ Thus, indistinguishable photon pairs have a slightly better concurrence than 
distinguishable ones, as expected from the visibility calculations provided in the previous section. 
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Figure 5: Concurrence as a function of average photon-pair number per pulse fi. Solid line, distinguishable pairs 
(Eq. |rj5|) ). dotted line, indistinguishable pairs (Eq. ([71)0. 

5 Correction terms by dark count 

From the results described in the previous sections, it is clear that visibility and concurrence can be made 
infinitely close to unity simply by decreasing /x. However, the dark counts of detectors limit those values in real 
experiments. Here we include the dark count in the theory developed in the previous sections. 

We assume that x photons are incident on a threshold detector with a collection efficiency a and a dark 
count probability d. When a and d are both much smaller than 1, the detection probability is approximately 
given by xa + d. With this approximation, we can calculate the coincidence rates in the previous sections by 
replacing xa with xa + d. 

In general, the two detectors have different collection efficiencies and dark count probabilities. We now 
denote their signal and idler values with subscripts s and i. 

For example, the coincidence rate of an HH measurement on indistinguishable pairs is given by 





{(x - y) 2 a s ai + (x - y)a s di + (x — y)aid s + d s di} 




(75) 



21 



The first term corresponds to the result shown by Eq. • The second and the third terms give the coincidence 
between counts caused by a photon and a dark count, and the fourth is the coincidence between dark counts. 
With a similar procedure, we can include the dark counts for all the above equations. In Appendix [Cj we 
summarize the equations that take account of the dark counts. We can optimize n to obtain the maximum 
visibility and concurrence using these equations. 

6 Time-bin entangled photon pairs 

In this section, we derive the equations for time-bin entangled photon pairs 4J. A time-bin entangled state is 
expressed with the following equation. 

|*) = -^(|l) a |l) i + |2) a |2)0 (76) 

Here, \k) z denotes a state in which there is a photon in the fcth time slot and in a mode z. We can generate 
a time-bin entangled photon pair by pumping a nonlinear optical medium with a coherent double pulse [H 
HI HH]. Therefore, the interaction Hamiltonian of the spontaneous parametric process that generates time-bin 
entanglement is exactly same as Eq. (j4]), if we replace the H and ^polarization modes with the first and second 
pulses, respectively. As a result, the pair number distribution of the indistinguishable time-bin entangled photon 
pairs is the same as that of the polarization entangled photon pairs given by Eq. (|19p. Similarly, the number 
distribution of the distinguishable pairs is Poissonian as given by Eq. ([25]) . 

In the measurement, each photon is passed through a 1-bit delayed Mach-Zehnder interferometer. If we 
denote the two output ports of the interferometer as a and b, the state \k) z can be transformed by the interfer- 
ometer as follows. 

\k)z -»i(|fc,o),-|*,6), + e M -|fc + l l o), + e w '|* + l,6),) (77) 

Here, \k, x) z indicates a state where there is a photon in the fcth time slot, the output x, and the mode z. With 
this equation, Eq. (|76p evolves into the following state by passing through the interferometers. 

I*) -> ^{|l,o) a |l,a) < -|l,o) a |l,6) < -|l > 6) a |l,o) i + |l,6) a |l,6> i 
+{ e W.+*i) + 1}|2, a) s \2, a), + { e *(»-+«0 - 1}|2, a) s \2, b) t 
+{ e <(*.+fl«) _ i}| 2 , fe) s |2, a) t + { e W'+ fl *> + 1}|2, 6>,|2, 6)< 

+e l ^+^ |3, a) s |3, a), + a ) s | 3) b) t + e <(«.+«0j 3) & ) s | 3 , a )i + e^ '+ e ^\3, b) s \3, 6),} , (78) 

Here, the terms that do not contribute to coincidences are not shown. When only ports a are connected to 
photon detectors, a constructive TPI occurs when 8 s + 9i = 2nn (n: integer). With a 100% collection efficiency, 
the coincidence rate per pair is 1/8. Therefore, with the same average pair number, the coincidence probability 
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of the time-bin entangled photon pairs is 1/4 of that of polarization entangled pairs. Here, note that the average 
pair number for the time-bin entangled photon pairs is defined as the number per two temporal modes. The 
decrease of the coincidence probability comes from the fact that the loss for each photon effectively increases 
by 3 dB because the photon is observed in the 2nd time slot with a 50% probability. Therefore, by replacing a 
with a/2 in the equations in section 3, we can obtain the formulas for time-bin entanglement. 

Let us describe an example. In a TPI measurement of a time-bin entanglement, the peak of the fringe 
is obtained when the detectors are both set at port a or b. So we denote the peak coincidence rate as R aa - 
Similarly, the coincidence rate at the bottom of the fringe is denoted as R a b- With the presence of dark counts, 
the TPI peak coincidence rate of the indistinguishable time-bin entanglement Raaindis is calculated as 



Raaindis — ^ ^ 



, ^ X+l 

x=0 \y=0 



1-1 



a\ x ~v 
2 



(79) 



which is 1/4 of RnHindis given by Eq. (|43|) . With an approximation 1 — (1 — a/2) x ~ xa/2 and with non- 
negligible dark counts, R pea kindis is obtained as 



OO I X 



R<i 



aindis 



x=0 Ky=0 



4 4 

The coincidence rates for other cases are summarized in Appendix! 



(LgL; 



(80) 
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Figure 6: Setup for time-bin entangled photon-pair measurement. 



7 Summary 

We described a theoretical analysis of the quality of entangled photon pairs generated by a spontaneous para- 
metric process. Assuming the use of threshold detectors, we derived simple formulas for the TPI visibilities for 
both distinguishable and indistinguishable polarization entangled photon pairs. The main results obtained in 
this paper are summarized below. 
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1. We presented the exact formulas for the coincidence rates of the peak and bottom of TPI fringes, taking 
account of the accidental coincidences caused by all the high-order multi-pair emissions. With these 
formulas, we can calculate the visibilities as a function of the average pair number and the photon collection 
efficiencies. In addition, these somewhat complex formulas can be reduced to simple functions of the 
average pair number when the collection efficiencies are small. These approximated formulas also include 
the effects of all the high-order multi-pair emissions. 

2. Using the derived equations, we compared the visibilities of the distinguishable and indistinguishable 
entangled pairs. The result suggests that the indistinguishable pairs provide improved visibility with the 
same average pair number, implying that we can construct quantum information systems with better 
performance using the indistinguishable entangled photon-pair sources. 

3. The full density matrix is obtained as a function of the average pair number and the collection efficiencies. 
Here, too, simpler formulas for the density matrix were presented. In addition, we derived simple formulas 
of concurrences as a function of the average pair number for distinguishable and indistinguishable pairs. 

4. We included the effect of dark counts in the theory, which enabled us to predict the visibility experimentally 
obtained with a specific average photon number. 

5. We modified the above formulas for time-bin entangled photon pairs. 

We believe that the contents of this paper will be useful for designing and evaluating photonic quantum 
information experiments with entangled photon-pair sources based on spontaneous parametric processes. 
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Appendices 

A Coincidence-to-accidental measurement of photon-pair state 

This appendix describes the coincidence-to-accidental ratio (CAR) measurements, which is used to evaluate 
correlated photon pairs (i.e. a source that generates a product state but not an entangled state) [2H1[27]. We 
derive the theoretical CAR for distinguishable and indistinguishable correlated pair sources. 
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Figure 7: (a) CAR measurement setup, (b) typical histogram of coincidence events. 



A typical CAR measurement setup is shown in Fig. [7| (a). A correlated photon pair source outputs signal- 
idler photon pairs, which exhibit a temporal correlation. Each photon is detected by a threshold detector. 
The detection signal from the detector for the signal photon is used as a start pulse for the TIA. The idler 
detection signal is delayed by an electrical delay circuit and then used as a stop pulse. We insert the delay 
so that the coincidence events caused by the photons generated at the same time do not fall within the TIA 
deadtime. Then, we can obtain a histogram of coincidence events as a function of time, as shown in Fig. [7] (b). 
If there is a temporal correlation between a signal and an idler photon, we observe a peak in the histogram in 
the temporal position that corresponds to the coincidence events caused by the photons generated at the same 
temporal position. We call the time slot where this peak is located a "matched slot". We also observe other 
side peaks, which are the coincidence events induced by the photons generated at different temporal positions. 
We call these time slots "unmatched slots" . 

We assume that the probability of obtaining x pairs is given by P(/i,x). Then, the count rate per pulse at 
each detector is given by 

oc 

Ci(/i) = ^P(/i^){l-(l-ar}. (81) 

X 

The coincidence rate at an unmatched slot R um is simply given by the coincidence probability of two uncorrelated 
events, which is expressed as 

Rum= (£p(M^){l-(l-af}) . (82) 
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On the other hand, the coincidence rate at the matched slot R m is given by 

oo 

i^ = X>(M,*){l-(l-a) x } 2 - (83) 

X 

The CAR C is given by 



C=|^. (84) 



A.l Distinguishable pairs 



When the source generates distinguishable correlated pairs, it is easy to show that the pair number distribution 
is Poissonian and given by Eq. (|23p as in the case of entangled pairs, using a similar procedure to that described 
in the previous section. We can obtain the CAR for distinguishable pairs simply by replacing P(/i, x) with 
P p (fi,x) in Eqs. (JUJ), (83]) and fMjl . 

When the collection efficiency is small, C can be approximated in a much simpler form. If a << 1, the 
detection probability of x photons {1 — (1 — a) x } can be approximated with xa. Then, using Eq. (|94p . R um 
can be simplified to 

/oo \ 2 

Rum — ^ Pp{n, x)xa} =i?a 2 . (85) 
Also, using Eq. ([95]) . R um is transformed to 

oo 

Rm ~ ^2P p (n,x)x 2 a 2 . = a 2 /x(/x+ 1) (86) 

X 

Thus, CAR for distinguishable correlated pairs C v is given by 

C p = l + -. (87) 
This result agrees with the result obtained with the simple Poissonian model introduced in [27] , 

A. 2 Indistinguishable pairs 

When the source generates indistinguishable correlated pairs, the pair number distribution is slightly different 
from that for the indistinguishable entangled pairs given by Eq. (|19p . Using the interaction Hamiltonian given 
by Eq. |T]) and a similar procedure to that described in section 2.1, the probability of generating x pairs Pth(x) 
is given by 



P th (x) = — = ,2.1 (88) 



Y 2x _ tanh 2 * xt 
TP ~ cosh 2 xt 

It is straightforward to show that Pth{x) can be rewritten with the average pair number \x as 
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Thus, the number distribution of indistinguishable correlated photon pairs is expressed as the thermal distri- 
bution. The CAR for the indistinguishable pairs is obtained by plugging P t h(n,x) into Eq. (|84p. 
Let us consider the case when a << 1. R um can be approximated as 

R S ide ^ PthiP; x)axj =Li 2 a 2 , (90) 

where we used Eq. . Similarly, using Eq. (pT7| found in Appendix [BJ we can calculate R m as 

oo 

i? m ~^P th ( M ,a;)a; 2 a 2 =a 2 (2 A i 2 + Ai ) (91) 

x—0 

Thus, CAR is obtained with the following equation. 

C th = 2 + - (92) 

With Eqs. ([87| and (|92|) . the CAR obtained with a indistinguishable correlated photon-pair source is larger by 
1 than that obtained with a distinguishable source. 

B Details of calculation in section [3] 

The pair number distributions Pi n dis(^,x), P p (p,x) and Pth(^,x), which are defined by Eqs. (TT9"|) . (f2"3")) and 
([59")) , respectively, satisfy the following equations. 

oo oo oo 

^Ppifl, x) = y^PindisiV, X) = },Pth(lJ>, X) = 1 (93) 
x—0 x—0 x—0 

oo oo oo 

^ a;-P P (M> x) = xP indls (fi, x) = xPthip, x ) = M (94) 

IE— 21—0 21—0 

OO 

^ x 2 P p ([i,x) — fi + p? (95) 

OO g 

2J x 2 P indis (p, x) = /i + -^ 2 (96) 

OO 

^VP tft (M,z) =^ + 2^ 2 (97) 



z=0 



The derivation of Eq. (|3T)|) is given in the following 

1 

2 X y\(x — y)\ v ai x\ 



R HHp ~ « 2 ^ K]_- i7 ^_ Tj (^ y ) 2 ^e- 

x=0 Lj/=0 



2 — « oo[ T a; 1 
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When x > 1, we can use the following relationship. 



A x-y 2 x ~ 1 x 2 X ~ 2 

^ y\(x -y-l)\~ Jx~iy. (x-2)V ( } 



Using the above, Eq. (|98|) can be transformed as 



x=2 
oo 



2 +a A,\2(^-l)! e * 4(,-2)! e 



2 2 ^1 (x-l)P 7 4 ^Wx-2)! 



e 



2 ^\ (x-iy. [ >j 2 \ ^ (x-l)lj 4 ^ l(x-2)! 
V 2 4 



(100) 



Also, the following relationships are used in the calculations in section [3] The first equation assumes that x 
is an integer that is larger than 1. 

x ^ 2 t 2 1 x ~ 2 

^ o (y~iy.(x-y-iy. = Jx~iy. (101) 

j2^-{x-yf = \x{l + 2x) (102) 
^-^ x + 1 6 

j/=0 

y ( * ~ v)v = x{x ' l) (loa) 

^ x+ 1 6 v ; 

C Coincidence rates of correlated and polarization entangled pho- 
ton pairs with dark counts included 

C.l CAR, distinguishable 

Rpeak — fJ-a s ai + ([ia s + d s )(fiai + di) (104) 

Rside — (M«s + d s ){pa l + di) (105) 

C p ~ 1 + - ^ — (106) 

(Lia s + d s )([iai + di) 



2S 



C.2 CAR, indistinguishable 



Rpeak — (J-a s ai + (i 2 a s ai + (fj,a s + d s )(fia l + dj) (107) 

Rside — (fJ,a s + d s )(fiai + di) (108) 

C 1 + ^ asQ:i (109) 

(fia s + d s )(/iOi + di) 



C.3 Entanglement, distinguishable 



-R/fffp ~ i/ita s a, + Q/-ta s + d s ^ + d *) ( 110 ) 

Rhv p ^ Q^ as + rfs ) Q^"* + 

- ^a s «. + + d s j Qm^j + ( 112 ) 

C.4 Entanglement, indistinguishable 

RHHindis ^ «s«i y + ^H 2 ' 2 

pi 2 aai jua a dj . /"Ms . , , 

RHVindis ^ T 1 2 1 2 1 

Rn+indis ^ a s "i I — + y I H ^ 1 2 

D Coincidence rates of time-bin entanglement with dark counts in- 
cluded 

Here, the average pair number /1 is denned per two time slots. In contrast, the average pair number per pulse 
was used in equations in [IS]. Therefore, the equations in [IS] provide the same results as the equations for 
distinguishable pairs shown below. 



D.l Distinguishable 



Raa P ^ ^na s ai + + d s ^j (y-^mi + di^j (116) 

Rab P ^ ^fia s + d s J \^Liai + d^J (117) 
R a +p ^ Yq M Q! sQ!i+ (^® s + d s j (-fiai+dij (118) 
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D.2 Indistinguishable 

Raaindis ^ a S «i ( ~§~ + g ) + "If" + "If" + ^ ( 119 ) 

/x 2 oai MOgrfj fia>jd s 

ftabindis — — 1 1 r «s«i K lM ) 

R a+indis - 0.0,(^ + ^+^ + ^+4* (121) 
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